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摘 要 复合 发 动机 预 冷 器 以 其 轻 质 紧凑 的 结构 且 高 效 的 传 热 性 能 逐渐 应 用 于 超 音 速 飞机 发 动机 中 , 研究 高 流速 高 温差 
极端 环境 下 微小 尺度 紧凑 结构 内 换 热 机 理 对 于 提高 发 动机 性 能 具有 重大 意义 . 本文 针对 Sabre 复合 发 动机 的 预 冷 器 ， 建 
立 三 维稳 态 可 压缩 流体 横 掠 又 排 管 束 的 强制 对 流 换 热 模型 ， 研 究 管 内 流体 速度 ,， 管 外 流体 速度 ， 入 射 角度 以 及 管 间距 对 于 
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管 外 空气 换 热 性 能 的 影响 机 制 . 结果 表明 , 通过 增 大 管内 制冷 介质 的 流量 , 缩小 管 间距 以 及 空气 垂直 入 射 等 几 种 手段 , 均 
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Abstract Micro-channel air pl is gradually used in supersonic aircraft engine, due to its 


t structure and high-efficiency in heat transfer. Thus, the numerical 


advantages of light weight, 2 
analysis of heat transfer ”> oler for hybrid airbreathing rocket engines at high flow rates and large 


temperature difference 


mathematical model is established to investigate the forced convectio 


great value to increase engine efficiency. In Fo paper, a three dimensional 


at transfer characteristics in 


the staggered tube bundle, in which the flow is steady, compressile, {orced convection. In particular, 
the effects of the flow rate, incident angle and the tube pitch Me analyzed. It is found that the heat 
transfer is enhanced when the coolant flow rate incre q the tube pitch decreases， As for the 
influence of incident angle, the results indicate that AS rtical incidence is better than other incident 


degrees. 


\S 
0 引 言 


超 音 速 飞 机 为 降低 飞机 自重 使 用 复合 式 预 冷 发 
动机 ， 通 过 预 冷 大 气 中 的 空气 作为 氧化 剂 ， 同 时 实 
现 吸 气 和 火箭 动力 提供 由 轻 质 大 功率 预 冷 器 作为 
发 动机 热力 循环 中 的 关键 组 件 ， 可 在 极 短 的 时 间 内 
从 来 流 空气 中 吸取 足够 的 热量 ， 达 到 深 冷 空气 的 目 
的 .因此 预 冷 器 的 功 重 比 与 紧凑 度 较 高 ， 常 采用 微 
小 通道 换 热 增强 其 换 热 性 能 。 国内 外 针对 复合 式 预 
冷 发 动机 的 研究 主要 包括 英国 的 Sabre -9， 日 本 
的 ATREX 上 5,6] 和 S-engine 三 种 发 动机 . 三 种 换 热 器 
虽然 结构 不 尽 相同 , 但 均 采 用 管 壳 式 , 紧凑 管束 的 换 
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热 方式 ， 沿 圆周 方向 布置 管束 ， 一 方面 使 换 热 器 结 
构 更 加 紧凑 ， 另 一 方面 通过 减 菏 管 壁 的 方式 增 大 换 
热 面积 ,减轻 换 热 锅 重 量 ， 

本 文 针对 Sabre 发 动机 预 冷 右 的 结构 进行 模拟 
验证 研究 。 由 于 预 冷 右 结 构 的 复杂 性 ， 对 于 该 预 冷 
器 换 热 性 能 的 研究 较 少 . 本 文 针 对 可 压缩 流体 横 护 
又 排 管束 的 强制 对 流 换 热 ， 建 立 三 维稳 态 模型 ， 研 
究 了 管内 的 流体 速度 ， 管 外 流体 的 速度 ， 入 射 角度 ， 
以 及 管 间距 对 管 外 至 气 换 热 性 能 的 影响 ， 对 预 冷 需 
的 结构 设计 具有 重大 指导 性 意义 ， 
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1 预 冷 器 几何 模型 简化 结构 中 , 用 两 个 半 管 和 一 个 完整 管子 , 体现 又 排 
结构 .对 于 两 个 半 管 的 横 截 面 , 选用 对 称 边 界 条 件 ， 
1.1 预 冷 器 结构 使 得 半 管 的 流动 状态 与 完整 管 相 似 .此 结构 一 方面 
预 冷 右 的 基本 结构 如 图 1 所 示 , 毛细 管 结构 如 。 尽 可 能 简化 了 模型 ， 隆 低 网 格 数量 ， 同 时 不 影响 计 

图 2 所 示 ， 算 结果 . 
此 外 ,如 图 3 所 示 , 对 管子 进行 编号 , 1~25 依 
六 组 入 组 次 沿 顺 时 针 分 布 . 分 析 可 知 , 25 号 管 在 A-A 截面 的 


预 准 片 ~ 流动 状态 ,与 1 号 管 在 B-B 截面 的 流动 状态 相似 . 

上 上 依次 类 推 , 24 号 管 在 A-A 截面 的 流动 状态 , 与 25 号 
管 在 B-B 截面 的 流动 状态 相似 . 归纳 可 知 ，A-A 截 
面 和 B-B 截面 之 间 的 不 同 管 可 以 等 效 看 作 1 根 完整 
的 管 ， 因 此 可 截取 1 号 管 到 2 号 管 的 部 分 作为 模型 
计算 基础 单元 , 如 图 4 所 示 . 将 图 4 所 示 单 元 的 管 
束 部 分 等 效 看 作 1 号 管 的 不 同 部 分 ， 即 认为 由 外 向 
内 的 管子 沽 %- 根 管 。 最 外 层 管子 的 出 口 参数 对 应 下 
ea 个 


Re 口 。 


图 1 预 冷 右 模块 


Fig. 1 Precooler module 


J 


XX、 图 3 毛细 管 分 布 情况 
毛细 管 SS Fig. 3 The distribution of the micro-channel 


图 2 预 冷 片 示意 图 
Fig. 2 The diagram of the precooling piece 


预 冷 片 呈 弧 形 分 布 ， 弯 曲 向 内 ， 按 一 定 角度 贺 
周 排列 . 预 冷 器 的 内 外 图 支管 分 布 贺 直径 分 别 为 340 
mm、430 mm。 弯曲 半径 取 190 mm, 沿 径 向 按 8 层 
布置 ，25 个 预 冷 片 围 成 一 个 基本 结构 .单个 预 冷 片 
毛细 管 数 量 设 定 为 80， 毛 细 管 外 径 1.2 mm, 壁 厚 1 
0.15 mm., 

1.2 模型 描述 与 简化 we 

由 于 管 外 为 又 排 管 束 ， 因 此 为 了 体现 又 排 结构 

特性 , 需 建 立 三 维 流动 模型 , 模拟 管内 外 换 热 结构 . 建立 几何 模型 如 图 5 所 示 , 管内 液 氮 由 内 圆 流 
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向 外 圆 ， 管 外 空气 沿 半径 方向 入 射 再 沿 轴线 方向 被 
抽 走 . 


图 5 简化 模型 结构 
Fig. 5 The simplified model 


E> 


李 晨 沛 等 : 复合 发 动机 预 冷 器 换 热 特 性 研究 


其 中 4n 和 hout 分 别 为 管内 外 的 壁面 面积 ， hs, 和 
hout 分 别 为 内 外 流体 的 换 热 系数 ，T，T! 分 别 为 内 
外 流体 接近 壁面 处 的 流体 温度 , TT,, 7 为 内 外 壁面 
温度 . 
2.3 沸腾 换 热 计算 公式 

在 考虑 沸腾 换 热 四 情况 下 , 对 流传 热 系数 计算 
公式 如 下 : 

Qb = Sanb 十 Qep (7) 

其 中 ， aw 为 平均 沸腾 传 热 系数 ，anb 为 泡 核 沸腾 传 
热 系 数 ，aceb 为 膜 核 沸 腾 传 热 系数 ，s 为 泡 核 沸腾 抑 
制 因子 。 


amb = C0 Mo pr (— lgpr) (8) 


其 中 ， 伦 0.33 
C= 90 W033/(m0.66 .区 ) (9) 
CN m=0.12= 0.21g|Rpli, (10) 


对 面相 粮 度 ， 民 p 二 0.3~0.4 um; ML 为 液体 的 分 
量 , g/mol; pr 为 对 比 压力 (液体 压力 与 该 流体 的 
临界 压力 之 比 ); q 为 热流 密度 ,W/m?. 


0.736 
和 oe HE be 2 其 中 ， Qt 为 液体 对 流传 热 膜 系数 ， Xstt 为 Marttinelli 
ap .Olpu) op Opw) _ 1 (1) 
Ot Oz Oy Oz LN 
O(pu) 0 人 ao 加 四 
所 +div(puU) = div(7 .gradu) 十 9u 一 区 (2) 
XX 至 LAL 
人 上 +div(pouU) = div(n .grado) + Sv — 村 (3) Ne 6 hoyf 疆 U0 
oO(pw) " ge pp ， 和 其 中 , pL 为 液 相 密度 ; om 为 汽 液 混合 物 的 均 相 平均 
页 de 一 dh sad 十 39 3 类 和 ”密度 ;hio 为 汽 液 两 相 全 部 为 液 相 时 的 换 热 系数 . 
此 外 ,本 文 考虑 换 热情 况 ,引入 能 量 方程 澳 ， 
a Olpuh) Opvh) Opwh) _ 5 二 TOS 10 0 Re (14) 
t Oz Oy Oz (5) 
—p:divU +div(A:gradT)+ Sh+w 其 中 ， 
式 中 ， U,V 和 WwW 分 别 为 流体 速度 矢量 U 在 三 个 坐标 Retp ReL F125 一 一 a 一 万 1.25 (15) 
方向 上 的 分 量 , p 为 流体 微 元 压力 , p 为 流体 密度 , 7) HL 


为 动力 莫 性 系数 , 5,, 5S, 和 Sv 为 广义 源 项 , 和 为 流 
体 导热 系数 ,5 为 耗 散 能 ，Sh 为 流体 内 热源 ， 本 文 
中 Sh=0, 
2.2 热平衡 

同时 考虑 管内 外 流体 流动 对 换 热 性 能 的 影响 ， 
保证 管内 外 的 热流 量 守 恒 有 


Q = Ainhin (Tr 本 Tw) Aput hout (Tf > Te,) (6) 


式 中 ，G 为 质量 流量 ,x 为 干 度 ,yi 为 饱和 液体 笑 
度 , D 为 管 径 . 
3 求解 方法 

流体 计算 分 管内 流动 与 管 外 流动 两 个 部 分 . 管 
外 空气 侧 流体 计算 采用 标准 k-e 满 流 模型 ， 侧 边 的 
边界 采用 绝热 边界 条 件 , 进口 采用 速度 边界 条 件 , 出 
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口 为 压力 边界 条 件 ， 与 壁面 连接 的 边界 保持 热流 量 
平衡 . 管内 液 氮 侧 流 体 计算 时 由 于 液 氮 流速 不 同 , 当 
管内 Re 数 小 于 2300 时 , 不 采用 消 流 模型 ,管内 Re 
数 大 于 2300 时 , 采用 标准 ke 湛 流 模型 ,两 侧 的 半 
管 面 均 采用 对 称 边界 条 件 ， 保证 整个 管子 结构 的 流 
动 对 称 性 , 进口 给 定 流量 和 温度 , 出 口 给 定 压力 . 同 
样 与 避免 连接 的 边界 保持 热流 量 平衡 。 为 保证 管子 
进出 口 的 连续 性 . 使 前 一 根 管 的 出 口 参数 与 下 一 根 
管 的 进口 参数 保证 连续 ， 
压力 速度 耦合 采用 SIMPLE 算法 , 动量 和 能 

均 采 用 二 阶 迎 风 差 分 格式 , 压力 、 密 度 、 体 积 力 、 动 
量 、 能 量 的 松弛 因子 分 别 为 0.3、1、1、0.7、1， 连续 
性 方程 、X 方向 分 速度 、Y 方向 分 速度 的 残 差 控制 
为 10, 能 量 方程 残 差 控 制 为 10-. 


4 模拟 结果 及 讨论 


4.1 液 氮 流 量 对 于 空气 侧 换 热 性 能 的 影响 


本 文 首先 对 沸腾 式 换 热 和 单 相 换 热 的 结果 做 SN 


了 比较 . 由 于 沸腾 式 换 热 的 管内 压力 在 临界 压 
下 ( 液 氮 临界 压力 为 Prit=3.3958 MPa)， 


虑 沸腾 式 换 热 时 ,设置 管 间距 $1=3.6 内 入 
口 压力 为 P=3 MPa， 入 口 温度 为 K， 流 速 


ua=0.24 m/s( 展 流 ). 空气 入 D473,=293 K， 
入 口 压力 Pir=0.1013 MPa “入 速 uair=10 m/s。 根 
据 沸 腾 式 换 热 计算 公式 计算 所 得 ， 空 气 出 口 温度 为 
Twir,out 二 273 KK， 空 气 侧 Nu=10.219. 

针对 单 相 换 热 ， 本 文 模拟 了 液 氮 流 速 在 0.24 
m/s、1.5 m/s、2 m/s、4 m/s 人 情况 下 ， 液 氮 入 口 温 
度 为 76 KK, 压力 为 5 MPa, 管 间距 $1=3.6 mm， 空 


气流 速 为 10 m/s, 空气 垂直 入 射 (空气 入 射 方 向 与 管 . 
内 制冷 介质 入 射 方向 垂直 )， 入 射 温度 为 293 K, xdfo' 


口 压力 为 1 个 大 气压 下 ， 空 气 侧 的 换 热情 况 XG 如 下 
图 6 所 示 . 

由 图 6 可 以 看 出 ,， 随 着 液 氮 流 量 的 升 高 ， 空 气 
侧 的 Nu 数 增高 , 换 热 性 能 增强 , 当 液 氮 流速 由 0.24 
m/s 增加 到 1 m/s 时 ，Nw 数 由 11.09 增加 到 11.31， 
换 热 性 能 增强 较 大 , 这 是 因为 在 0.24 m/s 的 情况 下 ， 
管内 液体 流动 状态 为 层 流 ， 换 热 性 能 较 差 ， 而 随 着 
液 氨 流 速 的 增加 , 到 1 m/s 以 后 已 转 为 漠 流 , 换 热 性 
能 增强 . 当 流速 增 大 到 1.5 m/s 后 , 换 热 性 能 的 增加 
量 减少 , 因此 , 建议 管内 液 氮 流 速 控 制 在 1.5 m/s. 

沸腾 换 热 与 单 相 换 热 性 能 比较 可 知 ， 沸 腾 换 热 
的 得 到 的 出 口 温度 为 Tui,out=273 K, 单 相 换 热 在 1.5 
m/s 下 得 到 的 出 口 温度 为 Tas,out==262.675 K, 因此 ， 
沸腾 换 热 的 换 热 性 能 比 超 临界 条 件 下 的 单 相 换 热 的 
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换 热 性 能 差 , 此 外 , 从 安全 性 能 的 方面 考虑 , 建议 采 
用 单 相 换 热 ， 


25 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
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液 氮 流量 变化 对 空气 侧 换 热 的 影响 

Fi ® e variation of Nusselt number and temperature 
A he tube with an increase in flow rate inside the tube 
时 人 射 角度 变化 对 于 空气 侧 换 热 性 能 的 影响 

为 验证 空气 侧 入 射 角度 对 于 空气 侧 换 热 的 影响 ， 

分 别 模拟 了 入 射 角度 为 60", 65°, 70°, 75°, 80°, 85°， 

90° 情况 下 , 液 氮 入 口 温 度 为 76 KK， 压力 为 5 MPa， 

液 氮 入 口 速度 为 0.24 m/s, 管 间 距 为 3.6 mm， 空气 

流速 为 10 m/s, 入 射 温度 为 293 K, 出 口 压力 为 1 个 

大 气压 下 ， Bi 如 图 7 所 示 ， 


1 1 1 1 1 


60 65 . 70 75 . 80 85 90 
9/°) 
图 7 空气 入 射 角度 变化 对 空气 侧 换 热 的 影响 


Fig. 7 The variation of Nusselt number outside the tube with 


an increase in air incident angle 


由 图 7 可 以 看 出 ， 随 着 空气 入 射 角度 从 60° 增 
加 到 90°, 空气 侧 的 Nw 数 从 9.6 增加 到 11.08, 换 热 
性 能 增强 ,因此 , 垂直 入 射 效果 最 佳 . 

4.3 空气 速度 变化 对 于 空气 侧 换 热 性 能 的 影响 

为 验证 空气 入 射 速度 对 于 空气 侧 换 热 的 影响 ， 
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分 别 模拟 了 空气 入 口 速 度 为 8 m/s、 10 m/s、 12 m/s、 
20 m/s 情况 下 ， 液 氮 入 口 温度 为 76 K， 压 力 为 5 
MPa, 液 氮 入 口 流速 为 1 m/s, 管 间距 分 别 3.6 mm、 
4.8 mm、6 mm, 空气 入 口 为 293 K, 出 口 压力 为 1 个 
大 气压 下 , 空气 的 换 热 情况 , 如 图 8 所 示 . 


8 10 12 14 16 18 20 
Ua/m:s—! 


图 8 不 同 管 间距 下 , 空气 流速 变化 对 空气 侧 Nu wi 
1 


Fig. 8 The variation of Nusselt number outside the t 
an increase in flow rate outside the tube at ao itch 


RE 
本 
热 效 果 增强 ， 不 同 管 间距 条 御 环 的 Nu 数 增加 量 幅 
度 几乎 相同 . 但 是 由 于 空气 流速 的 增加 ， 会 造成 空 
气流 量 增加 , 进而 使 空气 侧 的 出 口 温度 上 升 , 如 图 9 
所 示 . 
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图 9 不 同 管 间距 下 , 空气 流速 变化 对 空气 侧 出 口 温 度 的 影响 


Fig. 9 The variation of temperature outside the tube with an 


increase in flow rate outside the tube at different tube pitch 


4.4 管 间距 变化 对 于 空气 侧 换 热 性 能 的 影响 
针对 管 间距 变化 对 空气 侧 换 热 性 能 的 影响 ， 本 


文 模拟 了 空气 流量 为 10 m/s 时 , 换 热 管 轴 向 间隙 为 
3.6 mm、 4 mm、4.8 mm、6 mm 时 ， 空 气 入 口 为 常 
温 , 出 口 压 力 为 1 个 大 气压 , 液 氮 流速 0.24 m/s, 入 
口 压 力 5 MPa、 入 口 温度 76 K 时 空气 的 换 热 情况 ， 
如 图 10 所 示 。 
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< 图 10 管 间距 变化 对 空气 侧 换 热 的 影响 


Fig. 10 The variation of Nusselt number and temperature 


outside the tube with an increase in tube pitch 


由 图 10 可 以 看 出 , 随 管 间距 的 增加 , 空气 侧 的 
换 热 系数 降低 ， 换 热 性 能 减弱 。 由 图 8 同样 可 以 发 
现 管 间距 增 大 换 热 性 能 减弱 . 此 外 ， 从 图 9 可 以 看 
出 随 着 管 间距 的 增 如 ， 空 气 出 口 温度 随 着 空气 流速 
的 增加 的 提高 量 猎 浙 降 低 . 对 于 小 管 间距 而 言 ， 空 
气流 过 的 朗科 对 空 所 出品 湿度 的 用 较 大 ， 说 明 
凑 型 又 排 管 束 的 换 热 性 能 更 好 1， 减 小 管 间距 有 
助 妃 所 高 换 热 性 能 ， 


结 论 


针对 管 沉 式 换 热 器 的 结构 ， 我 们 针对 不 同 的 党 
侧 流体 流速 、 管 内 流体 流速 、 壳 侧 流 体 入 射 方向 、 管 
间距 等 对 于 空气 换 热 的 影响 做 出 了 模拟 研究 .研究 
得 出 以 下 结论 : 

1) 沸腾 式 换 热 的 换 热 性 能 比 超 临 界 条 件 下 的 单 
相 换 热 的 换 热 性 能 差 ， 因 此 建议 采用 超 临 界 条 件 下 
的 单 相 换 热 更 好 。 液 气流 量 的 增加 导致 空气 侧 的 换 
热 效 果 增强 ,但 是 增强 幅度 会 逐渐 减 小 ， 因 此 选择 
1.5 m/s 的 液 氮 流速 不 但 可 以 增强 空气 侧 的 换 热 , 同 
时 经 济 性 更 好 ; 

2) 针对 不 同 的 空气 入 射 角度 而 言 垂直 入 射 的 换 
热 效 果 较 佳 。 主要 原因 是 空气 流速 方向 与 换 热 管 液 
氮 流 动 方向 的 不 同 夹 角 可 以 分 解 为 垂直 入 射 方向 的 
换 热 和 顺 流 换 热 , 由 于 顺 流 换 热 效果 较 差 , 因此 , 70? 
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及 60° 下 的 换 热 效果 没有 90° 的 好 ; 

3) 空气 流速 的 增加 虽然 会 增强 空气 侧 的 换 热 效 
果 , 但 是 由 于 空气 流量 相应 增加 ， 并 没有 使 空气 的 
出 口 温 度 降低 ; 

4) 对 于 不 同 换 热 管 轴 向 间距 的 研究 , 3.6 mm 间 
距 下 的 换 热 效果 较 好 ， 即 紧凑 式 的 换 热 结构 更 有 利 
于 换 热 效 果 的 提升 . 
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